Composite boson dominance in many-fermion systems 



Fabrizio PalumbcQ 
INFN - Laboratori Nazionali di Frascati 
P. O. Box 13, I-OOO44 Frascati, ITALIA 
e-mail: fabrizio . palumboOlnf . inf n . it 

I recently proposed a method of bosonization based on the use of coherent states of fermion com- 
posites, whose validity was restricted to smooth structure functions. In the present paper I remove 
this limitation and derive results which hold for arbitrary interactions and structure functions. The 
method respects all symmetries and in particular fermion number conservation. It reproduces ex- 
actly the results of the pairing model of atomic nuclei and of the BCS model of superconductivity 
in the number conserving form of the quasi-chemical equilibrium theory. 



I. INTRODUCTION 

There are many finite and infinite fermion systems 
whose partition function at low energy is dominated by 
bosonic modes. This is always the case when, due to 
spontaneous breaking of a global symmetry, there are 
Goldstone bosons. The cfFcctive bosons of these systems 
can be charged (fermion number 2), like Cooper pairs in 
metal, nuclear and so called color superconductivity[l] or 
neutral (fermion number 0), like phonons in condensed 
matter and pions in hadronic physics. 

The structure of composite bosons can also be changed 
by varying temperature or control parameters: for in- 
stance strongly interacting fermionic atoms in magnetic 
traps can form at some temperature molecular pairs 
which however condense only at a lower temperature. 
Moreover at zero temperature molecular pairs can be 
transformed into Cooper pairs by tuning an external 
magnetic field which controls a Feshbach resonance. For 
such systems two features are most relevant: the role of 
molecular versus Cooper pairs and the contribution of 
noncondensed pairs to the ground state energy [2]. 

There is an immense literature concerning techniques 
to derive effective Hamiltonians for both types of effec- 
tive bosons. I will refer to all of them as bosonization 
methods. To put my work into perspective I will briefly 
review the most relevant step in the development of this 
subject, but with no attempt to completeness. 

The first approach to bosonization of which I am aware 
is due to Bogoliubov Q. After the BCS work on super- 
conductivity Bogoliubov reformulated their theory using 
the Frohlich Hamiltonian [^j of electrons interacting with 
lattice phonons and rederived all their results concerning 
ground state properties. He then mapped the Cooper 
pairs into effective bosons and found that their dynamics 
is described by the Hamiltonian of a superfluid system 
of elementary bosons he studied previously, whose ex- 
citations are called Anderson-Bogoliubov sound. These 
modes are related to a nonvanishing contribution of non- 
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condensed pairs to the ground state energy density. 

The presence of bosons external to the electron sys- 
tem, the lattice phonons, plays an essential role in Bo- 
goliubov's theory. But in atomic nuclei in which nucleons 
are assumed to interact via a nucleon-nucleon potential 
(without mesons), in gaseous systems of fermionic atoms 
and in the BCS model there are no external bosons. In 
particular in the BCS model the contribution of noncon- 
densed pairs to the ground state energy density must van- 
ish since the BCS solution for the ground state energy- 
density is exact in the thermodynamic limit and 
therefore collective excitations cannot coincide with the 
Anderson-Bogoliubov sound. 

The importance of superconductivity in atomic nu- 
clei was immediately understood by Bohr Mottelson and 
Pines and the method of BCS, which breaks fermion 
number conservation was adapted to atomic nuclei for 
which this symmetry is important. Immediately super- 
conductivity of infinite nuclear matter was investigated 
by many authors [8|, and much later it was suggested 
that because of the strong tensor force Cooper pairs 
in symmetric nuclear matter should have the deuteron 
quantum numbers Q . Subsequently a smooth transition 
from Bose-Einstein condensation of deuteron-like bound 
states at low densities and temperatures to BCS pairing 
at higher densities was considered [Toj . 

More recently the superconducting properties of ultra- 
small metallic grains have been investigated [ll| ■ In this 
context, like for atomic nuclei, it is important a theo- 
retical method which respects fermion number conserva- 
tion [12. 

There have been many attempts to reformulate the nu- 
clear Hamiltonian in terms of effective bosons. Beliaev 
and Zelewinski [isj made an organic theory of bosoniza- 
tion with some points common to the present method, 
but their expansion has problems of convergence and vi- 
olates fermion number conservation. Later on much in 
the spirit of BCS but in a phenomenological approach 
Arima and lachello [14] introduced two different compos- 
ite bosons, the s— (spin 0) and d— (spin 2) bosons. Their 
model, the Interacting Boson Model, proved extremely 
successful in reproducing low energy nuclear properties, 
but it has not been derived in a fully satisfactory way 
from a fundamental nuclear Hamiltonian. One of the 
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features of the model to be understood is why can one 
(mostly) restrict the boson space to the s- and d-bosons. 
A justification will emerge in the method I will present 
which, even though is found with a very specific fermion- 
fermion interaction, might have a more general validity. 
I must notice that the absence in the Interacting Bo- 
son Model of bosons with the deuteron quantum num- 
bers predicted in infinite nuclear matter might be due to 
the large size Q of the deuteron-like Cooper pairs (very 
roughly) estimated of about 17 fm, and the fact that in 
heavy atomic nuclei the valence shells of the protons are 
different from those of neutrons. 

In general there are many methods of bosonization. In 
1+1 dimensions there is a wealth of exact results [l^. 
These results have been extended fl6j to many dimen- 
sions: first by use of renormalization group transforma- 
tions the fermion space is reduced to a tiny shell around 
the Fermi surface, then the problem is essentially reduced 
to one dimension by considering only excitations normal 
to the Fermi surface. In this way many general proper- 
ties can be investigated, but in many cases the effective 
parameters introduced by renormalization group trans- 
formations are difficult to evaluate. 

Different approach in multidimensional problems are 
based on several recipes jl7| for mapping of a fermion 
model space into a boson space. Such methods respect 
fermion number conservation and in principle yield an 
exact solution to the problem, but in practice one has to 
perform a truncation in the fermion space related to a 
selection of degrees of freedom guided by physical insight 
and calculational convenience. One shortcoming of this 
procedure is the appearance of " intruders" , namely states 
which in spite of their low energy do not appear in the 
boson space generated by the mapping p^ . 

For many-body systems like fermionic gases in mag- 
netic traps the most common approaches are based on 
numerical simulations or the quasi-chemical equilibrium 
theory [l^,[13l, but a true bosonization to my knowledge 
has not been achieved. 

I developed an approach to bosonization [2l| in which 
an effective bosonic action is derived by evaluating the 
fermion partition function in a basis of coherent states of 
fermionic composites. Coherent states offer for compos- 
ites the same advantages they give for elementary bosons 
and fermions. 

After bosonization is achieved the fermion dynamics 
can be studied by functional or numerical methods 22]. 
Analytic calculations are also possible in this approach 
both in the path integral and Hamiltonian formalisms, at 
the price of an expansion in the inverse of the dimension 
of the fermion space. Such an expansion respects fermion 
number conservation and therefore can be used also for 
finite systems, and in fact its first application was to 
atomic nuclei. 

An important issue was left over in the original paper. 
It concerns a subtraction necessary to perform the expan- 
sion in the presence of a condensate. The expansion was 
done under the assumption that the structure function 



of the condensed bosons be almost constant, a restric- 
tion which can be fulfilled in some atomic nuclei but is 
not in most systems, including the BCS model. For this 
reason it was not possible to test the method on this the- 
ory. This limitation is removed in the present paper, and 
general results are derived for arbitrary fermion-fermion 
interactions and structure functions of the composites. 
Actually the formalism is more general, as it will be il- 
lustrated at the end of Section III, where the effective 
bosonic action is reported. But let me anticipate that 
it has been applied also to relativistic field theories, in- 
cluding gau ge theories, at zero [2^ and finite fermion 
number [2j|. In the nonrelativistic domain the present 
method of bosonization finds its potential applications in 
finite systems like atomic nuclei and small metallic grains 
and in systems with nonseparable interactions like BEC 
of fermion gases. 

The paper is organized in the following way. To make 
it reasonably selfcontained I included some material from 
Refs. 2l[, [231 . In Sections 2 and 3 I outline the approach 
and the derivation of the effective bosonic action [2l|. 
whose details can be found in Appendix|D] Section 4 con- 
tains the derivation, valid for arbitrary fermion-fermion 
interaction, of the effective boson Hamiltonian which in 
Section 5 is expanded in the inverse of the dimension 
of the fermion space. The leading order, given by Eq. 
ITfl has the form of the Bogoliubov model of a superfluid 
boson system. 

As a test of the method for finite systems I apply it 
in Section 6 to the pairing model of nucleons in a single 
j-level. The spectrum '25] is correctly reproduced in a 
form useful for an understanding and justification of the 
Interacting Boson Model of nuclear physics. 

Section 7 contains an application to infinite systems. 
All the results of the BCS model of superconductivity 
concerning the ground state are exactly reproduced in 
the number conserving form of the quasi-chemical equi- 
librium theory [l^. 

In Section 8 I summarize the results and conclude with 
an outlook. 

As already said in the present paper I do not deal with 
neutral composites of particle-hole type, but I will intro- 
duce some neutral composites constructed with charged 
composites. Particle-hole bosons will be discussed by an 
extension of the method in a future work. 

At last I want to mention that the present a ppr oach has 
been applied also to relativistic field theories [23| , includ- 
ing gauge theories. There are technical differences, due 
to ultraviolet divergencies in relativistic theories, which 
make necessary a different way of evaluating the effec- 
tive action of the composite bosons. But the essential 
strategy remains the same based on the use of coher- 
ent states of composites in the framework of the transfer 
matrix formalism which is close to the Hamiltonian for- 
malism of nonrelativistic theories. The method has been 
tested on models of fermions with large number of flavors 
and quartic interactions in 3-1-1 dimensions, exactly re- 
producing the gap equation for spontaneous breaking of a 
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discrete chiral symmetry and the mass of the efTective bo- 
son appearing in these models. Moreover the structure 
function of the condensed bosons has been determined 
for the first time. Its spatial part turns out, surprisingly 
enough, to be identical to that of the Cooper pairs of the 
BCS model. 



II. OUTLINE OF THE APPROACH 

Consider the partition function of a system of elemen- 
tary (non composite) bosons 



where 



Z = iv 



exp I -l(iJgjgjj^-;in) 



(1) 



where T is the temperature, /i the chemical potential and 
n the number operator. A sector of n particles can be 
selected by the constraint 



T— — In Z = n. 



(2) 



A functional form of Z can be found by performing the 
trace over coherent states [2^ 



la) = cxp (iif ^ |0) , 



(3) 



where K are the particles quantum numbers, their 
canonical creation operators and uk holomorphic vari- 
ables. Coherent states satisfy the basic or defining equa- 
tions 



aK\a) = aK\oi) 



(4) 



where ax are canonical destruction operators. In terms 
of these states we can write the identity in the Fock space 



d/i(a*, a)(a|a) "'"|Q;)(a| 



where 



dii{a*,a) = W 



K 



(5) 



(6) 



Using this resolution of the identity the trace in the par- 
tition function can be evaluated with the result (26j 

Z = y c?^(q!*, a) exp [— S'gjgjjj(a*, a)] . (7) 

In the above equation "Sgjgjj-^ is the action of the particles 

'5'elem = {-aj Vtat_i -f iJglgj^(at , at-i) 

t 

-Ha^at-i} , (8) 
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(/t + l ~ /*) : 



(9) 



and H^^Q^{a\ a) is the Hamiltonian in normal order. 

In a system of fermions whose low energy excitations 
are dominated by fermion composites I can restrict the 
trace to these composites. The restricted partition func- 
tion can be written 



Zc = tr 



Vc exp 



(10) 



where Vc is a projection operator in the subspace of the 
composites. By analogy to elementary bosons I assume 
for it the approximate expression 



2m 



where /3* , (3 are holomorfic variables and 
states of composites 



(11) 



coherent 



|exp|^X]/3,Btj; 



(12) 



The creation operators of composites with fermion num- 
ber 2 are 

= E (b\) 4. , (13) 

^ mi,m2 

while those of fermion number zero are 

= 4= E 4m s™. ■ (14) 

V"J — ^ / mim2 



The fit's are fermion creation operators, m and J, rep- 
resent all the fermion and boson quantum numbers and 
r2j is the index of nilpotency of the J-composite, which 
is defined as the largest integer such that 



(15) 



I will assume for simplicity the index of nilpotency equal 
to half the dimension of the fermion space for all the com- 
posites. The matrices $j are the structure functions 
of the composites, which must be determined by a varia- 
tional calculation. Only the solutions with high index of 
nilpotency are acceptable. 

I call the states |/3) coherent because they share with 
coherent states of elementary bosons the property of a 
fixed phase relation among the components with different 
number of composites. But the basic property of coherent 
states cannot be fulfilled. Indeed 



(16) 
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This is a consequence of the composites commutation 
relations, which are not canonical 



Bjc. 



(17) 



In states with a number of composites n << f2, the above 
equations can be approximately satisfied provided the 
structure functions are sufhciently smooth. Indeed in 
such a case the last term is of order n/ft. But in states 
with n ~ f2, it is not possible to satisfy them even with an 
absolute freedom about the form of the structure func- 
tions (which are instead determined by the dynamics). 
The best we can do [2l| is to satisfy them for states with 
n + k composites, for fixed n ^ Q and \k\ << fi. 

As stated in the Introduction, in the present work I 
will study only composites with fermion number 2. The 
properties of the operator V are reported in Appendix [Cl 
in the form derived in Ref. f23| . 

The trace can be exactly evaluated [2l| (see below) 
yielding a functional form of Zc 



Zc 



d(3*dp 

2711 



exp 



(-5eff(/3*,/3)). (18) 



The expression of the effective bosonic action S'gg is re- 
ported in the next Section. This result holds under the 
only physical assumption of boson dominance and the 
approximation adopted for V . 

In many-body physics it is often used the Hamiltonian 
formalism. The Hamiltonian of the effective bosons, Hb: 
cannot be read directly from the effective action, because 
S'gff (/?* , (3) does not have the form of an action of elemen- 
tary bosons. Indeed it contains anomalous time deriva- 
tive terms, anomalous couplings of the chemical potential 
and nonpolynomial interactions, which are all features of 
compositeness. Therefore it has been necessary to devise 
an appropriate procedure to derive Hb, which is given 
in terms of boson operators w , b, (not to be confused 
with the composite operators , B) satisfying canonical 
commutation relations, so that 



(19) 



= tr ( - — {Hb - lJ.Bn) 



/is is the boson chemical potential and n the boson num- 
ber operator. 

Hb has a closed form but, for a practical use, it is 
necessary to perform an expansion in inverse powers of 
the index of nilpotency fJ. 



III. THE EFFECTIVE BOSONIC ACTION 

The most general fermion-fermion interaction can be 
written as a sum of separable terms, so that the fermion 
Hamiltonian can be given the form 



Hp 



K 



(20) 



The one-body term includes the single-particle energy 
with matrix e, the fermion chemical potential and any 
single-particle interaction with external fields included in 
the matrix M 



ho 



Hf + M. 



(21) 



The matrices Fk are the form factors of the potential, 
normalized according to 



(22) 



In order to evaluate Zc I divide the inverse temperature 
in No intervals of size r 



T = 



Not 



(23) 



Then as shown in Appendix |D] the Euclidean effective 
action is 

SesiP* ' /^) = ^ E ^ti" { 7 [T^ + rR B^S/tB] 
+2RB^hB-Y,9K [{R-1)fIFk 



where 



'[RB'^fI) ^tT{RFKB) - RB^fI^RFkB 



h = ho-Y,9KFKFK 

K 



B 



(3-B^ 



R^ {I+B^BY 



(24) 



(25) 



(26) 



(27) 



Notice in the third line a trace inside the trace. The 
variables /3*,/3 are always understood at times t,t — 1 
respectively. Explicitly, for instance 



Rt = (iL + BlBt-i 



(28) 



SqQ has a global U{1) symmetry which implies boson 
conservation. 

The fermionic interactions with external fields appear 
in the bosonic terms which involve the matrix M (ap- 
pearing in h). 

The dynamical problem of the interacting (composite) 
bosons can be solved within the path integral formalism. 
Part of this problem is the determination of the structure 
matrices Bj which can be done by a variational calcula- 
tion. The present approach then shares two important 
features with variational methods: the restriction of the 
fermion space to a subspace, the space of the compos- 
ites, and the variational determination of the structure 
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functions. But unlike standard variational methods ex- 
cited states are treated at the same time and on the same 
footing as the ground state. 

Before showing a concrete way of application, I want to 
list some characteristic features of the present formalism: 

1) It can be used with interactions quartic as well 
quadratic in the fermionic fields. The latter ones include 
electromagnetic and phonon interactions in the nonrela- 
tivistic domain, and renormalizable relativistic field the- 
ories in 3-1-1 dimensions [I^IIJ]. In the application to 
these theories the main change is the replacement of the 
exponential of the Hamiltonian by the transfer matrix 

2) It can be extended to cases in which not all the 
fermions bosonize, like odd atomic nuclei, systems of 
fermions and composite bosons in quasichemical equilib- 
rium, and relativistic field theories at finite fermion den- 
sity. This latter case has already been studied for QCD 
at finite temperature and baryon density (2^ . and one 
can see that the procedure adopted can easily be applied 
to nonrelativistic systems as well 

3) It allows the treatment of different coexisting com- 
posite bosons, for instance Cooper and molecular pairs. 
In such a case the equations for the structure functions 
of the composites should have more solutions, possibly 
with the same quantum numbers 

4) It provides the structure functions of the composites, 
and therefore their effective coupling to other fields, like 
the electromagnetic field 

5) It allows the introduction of composites with quan- 
tum numbers different from those of the form factors of 
the potential. An example is given in Section VI: with 
a purely pairing interaction in atomic nuclei, one needs 
also bosons with angular momentum different from zero. 



IV. THE EFFECTIVE BOSON HAMILTONIAN 

As already said the derivation of the boson Hamilto- 
nian is not straightforward, because S^q differs in many 
respects from the action of elementary bosons shown in 
Eq.®. I notice that in S^q 

i) the time derivative term is not canonical 

ii) the coupling of the chemical potential (appearing in 
h) is also noncanonical 

iii) there are non polynomial interactions because of 
the i?-function. This function becomes singular, as it will 
become clear in the sequel, when the number of bosons 
is of order 17, reflecting the Pauli principle. 

Let us start by examining the features of compositeness 
when the number of bosons is much smaller than fi. Since 
the expectation value of /3* •/? is of the order of the number 
of bosons, in this case we can perform an expansion of 
logarithm and i?- function in inverse powers of fl. From 
the logarithm I get 

^ tr In [1 + TRB^\/tB] = ^tr (B^tB) 
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The first term can be made canonical by normalizing the 
boson form factors according to 



J,K- 



(30) 



The other terms are then of order Notice that the 

diagonal condition is only a matter of normalization, but 
the off diagonal one must be compatible with the dynam- 
ics, and a redefinition of the /3's can be necessary to get 
it. 

Expanding the i?-function I get the following couplings 
of the fermion chemical potential 



fj,F tr 



BlBt-i-\iBlB 



t-i) 



(31) 



Only the first term is canonical. However the anomalous 
couplings can be eliminated to order by a redefini- 
tion of the chemical potential [2l[ , so that in the case of a 
small number of bosons the Hamiltonian can be derived 
without difficulty. 

But when the number of bosons is of order f2, an ex- 
pansion of logarithm and i?-function can be performed 
only after an appropriate subtraction, as explained in 
Section V. I assume, and I will verify later, that after 
such a subtraction the anomalous time derivative terms 
be of order fl^^. Then Zc can be written in terms of an 
auxiliary Hamiltonian H' as a trace in a boson space 



tr exp 



1 

T' 



H' 



(32) 



H' is obtained from 5gg- by omitting the time deriva- 
tive term, and replacing the variables /3* , /? by corre- 
sponding creation-annihilation operators , b. These sat- 
isfy canonical commutation relations and should not be 
confused with the corresponding operators I3\B for the 
composites 

H'=: hrl2RB^hB-J29K [(^^l) F^Fk 

[ K 

+RB''Fl^tr{RFKB)~RBlFlRFKB | : . 



(33) 



The colons denote normal ordering and 



1 



B = —b-B'' 



R 



1 + ■Bb-B'' 



(34) 



tr [B^BB'fVtB] 



(29) 



From H' I will derive in the next Section the boson 
Hamiltonian Hb, Eq. ([55]) . I must notice that the case 
of a number of bosons much smaller than f2 cannot be 
retrieved from the above equations. The reason is that 
in the former case fip = 0{il), while — 0(f7") in the 
case n ~ ri. 
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V. THE EXPANSION 

A rather general way to perform the subtraction nec- 
essary for the fl^^ expansion is to write the operator R 
in the form 



(35) 



where 



K1K2 



r = 



1 + ^ E ^hK.BK^Bi^ 



K1K2 



(36) 



Notice that R can be expanded with respect to fj, while F, 
which however does not contain creation-annihilation op- 
erators, must be treated exactly. The parameters 'I'%^k2 
are related to the expectation values (pp[J>K2)- K is im- 
portant to observe that such expectation values do not 
break boson (and therefore fermion) number conserva- 
tion, but will in general break other symmetries, like 
rotational invariance in deformed atomic nuclei. Their 
determination allows therefore to study thermodynamic 
or quantum phase transitions like breaking of rotational 
symmetry by axial or triaxial shapes in atomic nuclei [131 . 

For the sake of simplicity in [2l[ I restricted myself 
to cases in which the structure functions are almost con- 
stant, which justifies a subtraction independent of them. 
The formalism was then tested in the case of the so called 
pairing model [1^. The ground state energy was ex- 
actly reproduced, but the spectrum of excitations was 
not studied because the boson Hamiltonian contains cou- 
plings of all the bosons among themselves and was not 
solved. As we will see such couplings are an artifact due 
to the inadequacy of that subtraction. 



A. Subtraction in the presence of an s-condensate 

In the present paper I remove the restriction that the 
form factors should be almost constant, and consider the 
case in which a condensation occurs in a single quantum 
mode, called the s mode, whose quantum numbers will 
be denoted by "zero" 



' K 



1K2 — ^Ki,0 5k2,0 



(37) 



By condensation of the s-boson I understand that the 
occupation number of this mode is of order fi. The terms 
appearing in t) can then be classified according to 

{h%-r^)BlB^+ J2 VkMBk.b],^ ~ 
blBoY^bKBl ^ Vn, (38) 



provided the sum of the occupation numbers of noncon- 
densed modes is much smaller than fl and the structure 
matrices are of order il". Under these conditions one 
subtraction is sufficient, otherwise more subtractions are 
needed. Expansion of time derivative terms in 5gg- ac- 
cording to this classification, neglecting contributions of 
gives 



1 

— tr 
2r 



In (1 + BjBt) - In (l + BlBt-i 



2n 

K1K2 



where 



CK1K2 — T Bk^B 



t 

Ki- 



rn 



(40) 



In the derivation of the above equation I assumed 
tr Cko = 0, for K ^ and I disregarded the fluctuations 
of the product PqPq, whose contribution is of order fi"^, 
so that terms of the form 



^^tr(Co/^i C0K2 )/3o/3o^t (/3ifi /3k2 ) 



(41) 



are total time derivatives and do not contribute to the 
action. Then replacing PqPo by n, the number of bosons, 
the temporal terms become canonical if I impose the nor- 
malizations 



1, 



2n 



:tr CkiK2 ^ -^CkiqCqK2 



SK1K2, 



(42) 



It is worth while noticing that the ground state wave 
function of free fermions is included in the class of states 
which are assumed to dominate the partition function, so 
that the variational evaluation of the structure function 
will tell whether an actual condensation will occur or not. 
Indeed if I choose 



(43) 



where 9 is the step function and mp are the fermion 
quantum numbers at the Fermi surface, for large the 
normalization gives 



1 



(44) 



and all the coherent states are dominated by the term 
with Up fermions 



(/3o^"" n 



(45) 



7 



B. n ^ expansion in the presence of an s 
condensate 

To proceed with the expansion of the boson Hamilto- 
nian it is necessary to know how the couphng constants 
gK scale with fl. For infinite systems, since $7 — > oo, to 
get finite energies we must require gK ~ Such a 

behavior is also acceptable for many finite systems and I 
assume it in the following. As a consequence the fermion 
chemical potential /.ip is of order H,'^. 

At this point it is convenient to introduce the notations 



no = 05 



bUo. 



(46) 



Neglecting contributions of order ^ I get 



^Vki,K2 ^Ki^k^o&o j + H.c 



(47) 



where Eq is a c- number and H.c. represents the Her- 
mitian conjugate of the operator in the round brackets. 
The coefficients £ and V are functions of the operator no, 
which explains the presence of normal ordering. 

In many cases the double sum over Ki , K2 reduces to 
a single sum, K2 resulting conjugate to Ki according to 
a one-to-one correspondence K2 — Ki 



A'#0 

H.C. 



(48) 



From their general expression, reported in Appendix [El 
we see that Ec ~ ^ while the operators £k,Vk take 
values of order . 

H' has the form of the Bogoliubov model of super- 
fluidity, with an important qualification to be discussed 
below. Because of the absence of terms involving three or 
four operators of bosons out of the condensate, H' can be 
approximately (see below) diagonalized introducing the 
phonon operators 



Ak = -^bKbl, = A=b^Kbo 



In terms of these operators 
H' ^Ec + SEc 



(49) 



(50) 



where 



SE, 



c 



£k) 



Ek = \IS 



'K 



(51) 



The average number of noncondensed bosons in the 
ground state is [11] 



(E^^Kb 



KI 



2 ^ 



Ek 



1 



(52) 



Approximating the commutation relations of the phonon 
operators by canonical ones introduces errors of order 
51^^ only if this expectation value is much smaller than 
fi. This condition is also necessary for the classifica- 
tion of Eg. ([55)) to hold. Now even though the operators 
£k and Vk take values of order 51°, so that also the 
phonon energies Ek are of this order, it is well possible 
that {^K^o^\cbK) — 0{il). In such a case a unique 
subtraction is not sufficient and farther subtractions are 
necessary. 

Now the qualification mentioned above. As reminded 
in the Introduction in his reformulation of the theory 
of superconductivity starting from the Frohlich Hamilto- 
nian of electrons interacting with lattice phonons, Bo- 
goliubov found an effective Hamiltonian essentially equal 
to that he studied previously for superfluid bosonic sys- 
tems as far as the bosonic excitations are concerned. The 
coefficients £k and Vk of this Hamiltonian are such [2^ 
that the excitation spectrum is 



E^ 



\K\, \K\ 



0. 



(53) 



As a consequence {J2K^o^^KbK),^Ec = 0{fl). Such a 
result cannot hold for the BCS model, because the ex- 
pression of the ground state energy density derived by 
BCS is exact Q. 



C. Determination of subtraction parameter and 
fermion chemical potential 

In the determination of the subtraction parameters I 
meet with a subtlety. H' commutes with the boson num- 
ber operator, so I can select sectors with a given number 
of bosons. But I am not guaranteed that these bosons 
carry fermion number 2, because of the noncanonical cou- 
pling of the chemical potential. This fundamental prop- 
erty can be enforced just exploiting the parameters in- 
troduced by the subtractions. Indeed, denoting by E'Q(n) 
the lowest eigenvalue of H' in the sector of n bosons, I 
require that E'f^{n) be the lowest eigenvalue for n = ^np 



d_ 

dn 



1 



E'oin) — 0, for n — —np 



(54) 



The above equation, together with the condition ^ on 
the fermion number, determines one of the parameters 
r and the fermion chemical potential fip as functions of 
the number of bosons, r — r{n),]lp = Jlp^n), ensuring 
that these bosons carry fermion number 2. The boson 
Hamiltonian in the sector of n bosons is finally 



HB{n) = H'{r,fip) + 2fipn. 



(55) 
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It depends on n explicitly and through the dependence 
on n of f,JI p. Therefore also the matrices Bj, the form 
factors of the bosons, will depend on n, namely on the 
number of fermions. 

Notice that H' provides a mapping of the fermion in- 
teractions with external fields 



c^Mc - 



1 



-tr 



(56) 



The ground state energy of the auxiliary Hamiltonian in 
a state of n bosons is 

E'a{n) = Ec +SEc + (n, o| : £o(no)^o : \n,o). (57) 

Taking normal ordering into account 

{n, o\ : £o{fio) no : \n, o) = £o("-o) ^^o (58) 

where 

£o[no) = £o("-o) + 2no5o£'oo^oooo- (59) 

TT-o is the number of condensed bosons in a state of n 
bosons and the quantities DkiK2... defined in Ap- 
pendix B. The constraints ^ and ([51)1 for subtraction 
parameter and chemical potential are 



dE'o{7i) 
dE'oin) 



-2|,--trCo„ 



dfip 



("-0 



dSEc 
dn 



^tr[rCo„] 



= 0. 



(60) 



I have two equations and three unknowns: /ii?, t"^, '^o- 
But I remind that the expansion holds under the 
assumption 



(61) 



otherwise further subtractions are necessary. I assume, 
and I will discuss later when this assumption holds true 



d 



5Ec^O{n''), —5Ec^0{n-^). (62) 
on 



Then the first equation with the normalization condi- 
tion (gj) gives 

r2=n + 0(f]°), (63) 
and using this result in the second equation I get 



1 

2f]' 



rtr 



O50-D00 = o(o-i) 



(64) 



which determines ^Lp. 



I can finally write the ground state eigenvalue of Hb 



Eo{n) = E'oin) + 2nJIp = <{ ^tr 



TBo eBl 



-ngaDl^^n + 0{n°). (65) 

To evaluate the excitation energies in a state with / 
phonons I need the expectation value 

(n, /I : £0(^0) "0 : I", /) = ^0(^0 - /)(n-o - /) 

— d 

^ £a{no) - - — [£o(rio) f^o] /• (66) 
ono 

According to Eq. lpn]) the coefficient of / is of order fi"^, 
so that to order fi^ this term does not depend on f. 

The excitation energies of states containing / phonons 
are therefore 



EkiK2...Kj — Ek^ 



E 



K2 



E 



(67) 



VI. FINITE SYSTEMS: THE PAIRING MODEL 
FOR NUCLEONS IN A SINGLE j-SHELL 

I consider a system of nucleons in a single j-shell, in 
which case the only fermion quantum number is the third 
component of angular momentum m and the index of 
nilpotency ft is equal to j -I- ^. Composite bosons are la- 
belled by the angular momentum L and its third compo- 
nent M: K = {L,M). The form factors of the potential 
are proportional to Clebsh-Gordan coefficients. With the 
normalization of Eq. they are 



{Elm) 



LM) nil i^^^2 



'^{jmijm2\LM). 



(68) 



It is not necessary to solve the variational equations for 
the boson form factors, since because of rotational in- 
variance they also must be proportional to the Clebsh- 
Gordan coefficients 



B 



LM 



M- 



(69) 



As a further simplification I restrict myself to the pair- 
ing model, namely I assume the single-particle energy to 
vanish and the interaction to be a pure pairing poten- 
tial {gx — for K ^ ). I will evaluate the excitation 
energies to order il^ but the ground state energy only 
to order fl since it was already evaluated to order in 

I can immediately find the normalizations 



1 



1 



ao = 



a.L = 



1 



(70) 



and evaluate the coefficient oi H' . I find that L — L,M = 
—M, and 

Ec = 2hniy^ + ngo v 
Eq = [2/l-5o(^^ + 2 + 2^/](l-J/) 
Elm = [2/(,(l-2po)+4r!.gopo(l-po)] 
Vlm = (-l)*^2[-2/i + r!5o(l-2po)] (71) 
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where po = nn/fl. Eqs. ([63)) and ((64|) give 



r2 



With these values 



-aoin-n) + o{n-'). 



(72) 



Vlm = 0(r!-i), 5lm = r!,go + 0(f^"'). (73) 

Since the number of L- modes is of order fl^, namely 
6Ec — 0(ri"), the effective boson Hamiltonian in the 
sector of 2n fcrmions is 



gon 



^90 Y] b[j^jbLM + o{n- 



I / 
LM 



(74) 



This is the correct result |25| in the required approxima- 
tion. It is to be noted that noncondensed bosons have 
all the same energy Elm — ^go, so that the spectrum 
depends only on their total number, and not on their dis- 
tribution in different L-states. This property makes con- 
sistent the truncation of the boson space to a few modes, 
and as a consequence we can restrict the boson model 
space to only one of them, the d-boson in the Interacting 
Boson Model. 

The decoupling of different noncondensed bosons is 
achieved by means of the subtraction of Eq. ((37| . With a 
subtraction independent of the structure functions all the 
modes are coupled, so that the above properties might 
have been found only after solution of the resulting bo- 
son Hamiltonian. 



VII. INFINITE SYSTEMS: 
SUPERCONDUCTORS AND THE BCS MODEL 

I consider an infinite system of fcrmions whose quan- 
tum numbers arc spin and momentum, m = {s,p), and I 
parametrize the potential form factors according to 



(75) 



For simplicity I restrict myself to bosons of spin zero 
and momentum K and parametrize their form factors 
according to 



K 



Sl,S2"pi+P2,K Vpi 



K^-hK{K). (76) 



Since form factor and structure function are fully anti- 
symmetric 



f,{K) - /_,(i^), cj>,{K) = 0_,(if). 



(77) 



The structure functions are subject to the normalizations 
of Eqs.gSl) 



1 



(78) 



where 



(79) 



For infinite systems neglecting terms of order ^ in the 
evaluation of the ground state energy density does not 
produce any error in thermodynamic limit. So I will ne- 
glect such terms. Therefore the ground state energy is 



£^0 



3Er 



^tr[rSoeBt 



90 



ngoD't 



n 

n' 



(80) 



The expressions of 5Ec,£k and Vk can be easily eval- 
uated inserting the definitions of potential form factors 
and boson structure functions in the equations reported 
in Appendix B. 



A. BCS model 

For a general potential SEc ^ 0(ri"), and as already 
said consistency of the expansion requires further 
subtractions. But a great simplification must arise for a 
pure pairing interaction, gx = for K ^ 0, because then 
as reminded previously the contribution of 6Ec to the 
energy per particle must vanish in the thermodynamic 
limit, or in other words SEc must be of order ff. There- 
fore I assume SEc = 0{^^) and I will verify a posteriori 
that this is true. I also assume 0p to be real and I mini- 
mize the ground state energy imposing the normalization 
constraint by the Lagrange multiplier A 



d 



90 



^ ^ ^p^pfp 



I thus get the gap equation 



where 



(cp - A) yi^^p = - A (1 - (l>l)fp 



A = goV^^^p 4'p fp 



0. (81) 



(82) 



(83) 



is the gap function. These are exactly the results of 
the quasi- chemical equilibrium theory, Eq. IV (1.17) of 
ref. fed /. In this connection I remind that this theory, 
unlike the BCS theory, is fermion number conserving like 
the present approach. 

Following [29|] I then find A — epp,pp being the Fermi 
momentum of the system. Next I evaluate the coefficient 
of the quartic boson term 



p~K 



[2(e 
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P (1 - '^(t>p) ^ fp <Pp-K 

x<^p_i^(i^)</)p_i;,(-if)} + 0(r!-i). (84) 

In the above equation there appears the fermion chemical 
potential Jip which is determined by Eq. ([64]) . Using the 
gap equation I find Jip — A. Then using again the gap 
equation in the expression of I find that Vk = 0{il~^) 
and therefore SEc is at most 0{^lP). For a pure pairing 
interaction the contribution to the ground state energy 
per particle of noncondensed pairs vanishes, in agreement 
with 0/, and my assumption about 5Ec is justified. 

To get analytic expressions I assume to be the free 
fermion energy and the form factor of the potential to 
have the schematic form of BCS 



2m 



p\ fp^9{2mLo-\p^-pl\). 



(85) 



6 is the step function and lo a cutoff energy (identified 
with the Debye cutoff energy of lattice waves in metal 
superconductivity ). 

The solution of the gap equation is [l^ 



1 



1 



where 



(86) 



(87) 



VIII. SUMMARY AND OUTLOOK 

I extended the original formalism for boson domi- 
nance based on coherent composite states. The extension 
amounts to perform a subtraction proportional to the 
structure function of the condensed boson. The results 
obtained have a validity restricted only by the assump- 
tion that condensation occurs in a singlet state and that 
the coupling constants scale according to gx ^ ri^^. In 
particular, since fermion number is conserved, they hold 
for finite systems. Therefore the present bosonization 
method can find applications not only in the presence of 
important nonpairing interactions, but also to account 
for finitcness effects in atomic nuclei and small metallic 
grains. 

The general expression of the bosonic Hamiltonian con- 
tains a classical part SEc which couples the form factors 
of all the bosons. The theory becomes much simpler for 
pure pairing interactions, because SEc is negligible, and 
I tested it on two paradigmatic models. 

As an example of finite systems I considered the pairing 
model of nucleons in a single j-shell and I reproduced 
its bosonic spectrum in a form useful to understand and 
justify the Interacting Boson Model. Indeed I find that 
a boson space made of the s- and d-bosons is sufficient 
to reproduce the spectrum, which is not changed by the 
addition of other bosons. For the fermionic part of the 



spectrum one must use the technique developed in [24| . 
Needless to say the bosonization of the pairing model 
has been obtained in many ways, but I emphasize that 
the present formulation holds in a framework valid for 
arbitrary fermion- fermion interactions. 

I then applied this method to the BCS model of super- 
conductivity reproducing exactly its ground state proper- 
ties in the form of the quasi-chemical equilibrium theory. 

Two important issues are left for future work. The 
first one concerns the determination of structure func- 
tions for a general interaction and noncondensed bosons. 
I keep in mind in this connection the possibility of differ- 
ent structure functions with the same quantum number 
K . They might for instance be associated with intruders 
or coexisting molecular and Cooper pairs. 

The second one concerns phonons, namely neutral 
bosons describing polarization effects |30|. They must 
necessarily be included in many cases |3l| and certainly in 
the presence of particle-hole terms in the fermion- fermion 
interaction. 



APPENDIX A: BASIC FORMULAE IN BEREZIN 
INTEGRALS 

The definition of the Berezin integral for a single Grass- 
mann variable is 



^7(07 + b) = a, 



(Al) 



the generalization to many variables being obvious. For 
a change of variables 



in a multiple integral we have 



(A2) 



/n(^7.)/(7)= (detl^) ' j\{{di,)f{i). (A3) 

Notice the appearance of the inverse of the jacobian, con- 
trary to the case of ordinary variables. 

Gaussian integrals can be evaluated exactly, like for 
ordinary variables. There are two types of such integrals 

/ X{{d^ld-fn) exp ^ 7,* M.j7, = det M (A4) 
jY[id-fh)expY,ll^^^nj = Pf^ (A5) 



where Pf A is called the pfaffian [S^ of A. The follow- 
ing algebraic identity holds 



(Pf^) =det^. 



(A6) 
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APPENDIX B: INNER PRODUCTS OF 
COMPOSITE STATES 

Let us consider the case of only one composite. To 
evaluate the inner product of coherent states I use the 
identity operator in the fermion Fock space 



(Bl) 



where the 7*, 7 are Grassmann variables and I7) coherent 
states :32,1 



I7) = exp(-7 €■!■)). 



(B2) 



I then have 



^7*^7 exp(-7*7)(/3i|7) (7I/3). 



(B3) 

The matrix element (/3i|7) can be evaluated using the 
defining property of coherent states 



with the result 



{fii\l) = exp 



(B4) 



(B5) 



Therefore (/3i|/3) becomes 

(/3i|/3) = j dj*djE{j*,j,(3l,P), (B6) 
where the function E is 

E{r,j,(3*,P) = exp (^-7*7+^/3* 7 5 7 



+ 



2V^ 



By the change of variables 

p 

the integral is factorized according to 

(/3i|/3)= j d7'exp(^^7'/3i3ty^ x 



(B7) 



(B8) 



d-y exp 

The factors are of the form IA5[ so that finally 
1 



.(B9) 



1 



det ( — /3Bt 



det 



PIB 



det 



1 



(BIO) 



APPENDIX C: THE OPERATOR V 

In this Appendix I show that the operator V of EgfTTl 
approximates Vc to leading order in an expansion in the 
inverse of the index of nilpotency. 

For the sake of simplicity I consider the case of a unique 
composite. Then 



(CI) 



n=0 



where 



i.„- (0|B"| (#)"|0). (C2) 

I must then prove that 

(0|B™P(Bt)»|o) ^ (0|S"|(Bt)"|0) = <5,„,,j.„, . (C3) 
These equations are generated by the following ones 

{f3'\V\(3") ~ (/3'|/5") (C4) 



taking derivatives with respect to (}'* , (3" and setting 
these variables equal to zero. To simplify the formulae I 
adopt a slightly different definition of composites 



(C5) 



The right and left hand sides of Eq. IC4I are 

{f3'\V\P") = /^exp[-£(/3*,/3,r,/?")] 
iP'lP'') = exp[Trln(l+/3'*/3"BSt)] (C6) 

where 

£{f3* , /?, /?'* , /3") = Tr [ln(l + I3*I3BB'' ) 

-\n{l + (3'* (3 BB^) -\n{l + f3* 13" BB'')] . (C7) 

I evaluate the integral by the saddle point method. The 
saddle point equations are 



(/3 - /3")Tr- 







(1 + /?*/3BBt)(i + p*p"BB'f) 
, BB'^ 
{13* - r )Tr ^ ^*^5^t)(i + f3'*i3BB^) = ^'^^^ 

with solutions 

= ^p'* . (C9) 

At the saddle point 

£(^*,P,f3'*,f3") = -Tihi{l + f3'*P"BB^). (ClO) 
Moreover 



-Qp^\-^=t3",r=fi" "'^'"(i + /3'*/3"BBt)2 



: 

(Cll) 
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In conclusion 



where 



Tr- 



(C12) 



(1 +/5'*/5"SBt)2_ 

The desired result follows if we assume 

Tr {B^bY - ri-"+i . (C13) 

It is then easy to prove also the idempotency property of 
projectors 

V^V^. (C14) 



APPENDIX D: DERIVATION OF THE 
EFFECTIVE ACTION 

For the following manipulations we need the Hamilto- 
nian in antinormal form 

H = itr(/.+/io)-c/i^ct-^g^i cf^xci ct^ct (Dl) 

where the upper script T means "transposed" and h was 
given in Eq. (|25p . Now we must evaluate the matrix el- 
ement exp(— TiJ)|/3f_i). To this end we expand to 
first order in r (which does not give any error in the final 
r ^ limit) and insert the operator V between annihi- 
lation and creation operators 

{Pt\ exp(~TiJ)|A_i) = exp (^-^tiih + ho)T^ {/3t\V 
-ch^rVc^Y.9kr \ cFkcV ^ Flc^\l3t-i). (D2) 

k 

Using the identity in the fermion Fock space we find 



X exp ^— -tr(/i + hQ)T — j*hTj 
X exp I ^ gKT ^7 Fk 7 ^1*fIi* 

\ K 



(D3) 



where the function £'(7*, 7, /?*, /?) is defined in (|B7p . By 
means of the Hubbard-Stratonovich transformation we 
can make the exponents quadratic in the Grassmann vari- 
ables and evaluate the Berezin integral 

(/3t| exp(-ri7)|/3t„i) ^ / ]^ (ia|f dajf exp(-a* - a) 

K 

jitrln 1+ + 



X exp 



X detp exp ( — -tr(/i + /io)t 



(D4) 



p=l + hT. 



(D5) 



Performing the integral over the auxiliary fields a^* , 

we get 

(/3t| exp(-Ti7)|/3t_i) = j Y\da]^daK eyi-p{-a* ■ a) 



X exp < — tr In 



1 



X det p exp ( — 2^''^^^ + ho)T ) , 



(D6) 



The functional form of the composites partition function 
is 



Zc 



dp* dp 



2m 



exp(-5g5(/3*,/3)) (D7) 



where is given in Eq. 



APPENDIX E: GENERAL FORM OF THE 
COEFFICIENTS IN THE AUXILIARY 
HAMILTONIAN 

I will use the definitions 

TK1K2...K21-1K21 = -^^^[CkiK2---Ck21-3K21-2 



DKiK2...K2,-iK2 



xTBK2,^,hBl,^ 

— tr [CkiK2---Ck21-3K21- 

xTBk2,.,fI^, 



(El) 



for I an arbitrary integer. Let me emphasize that only 
the last factor in T involves the single particle kinetic 
energy h, and only the last factor in D involves the form 
factors of the potential. In the following I neglect the 
coupling to external fields and I assume that if only 2 
indices if 1 , K2 are different from zero, they must be equal 
in the £- terms. 

The classical energy is 



K 

The s-boson energy is given by 

^0 = Tqo — ly Toooo + ^fjo 

X-Doo-Doooo] + 
1 



CoqF]^Fk 



(E2) 



|2?oo|^ + 2(^-. 



K 



^5^2^*'^ ^B^Fk^FkBI 



(E3) 
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The terms V are separated in their contributions from 
tlic kinetic and potential terms in the fermion Hamilto- 
nian 



= E 



kin 



E 



pot 

KtK2 



kin 



Vki,K2 = ^KuK2+^K^.,K2 



.pot 



whose expressions are 
kin 



E 



c.pot 

■-K1K2 



— TK1K2 ^ '^TKiK2oa 



[TK1OOK2 
—>^{ToK2KiOQa + ^^100X2 00] ^ 



Tr 



OK2K1O 



2 9oD, 



00 



no 



-|;kin 



(E4) V 



pot 

K1K2 



DK2K1OO + D0K1K2O — ^ {D0K1K2000 



no 

2n 



^2 [TKiOKia — i^IifiOKaOOo] 



K2mK . 



2 <^ ngoDo 



DK20K10 — ^ {DQK10K200 



+DkiOK200o)^~^ + ^2 



'Jl OKiOK^K20aK 



(E5) 
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